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Abstract 

We provide an equation of state for hot nuclear matter in /3-equilibrium by applying 
CNJ \ a momentum-dependent effective interaction. We focus on the study of the equation of 

state of high-density and high-temperature nuclear matter, containing leptons (electrons and 
muons) under the chemical equilibrium condition in which neutrinos have left the system. 
The conditions of charge neutrality and equilibrium under /3-decay process lead first to the 
evaluation of proton and lepton fractions and afterwards of internal energy, free energy, 
pressure and in total to the equation of state of hot nuclear matter. Thermal effects on the 
properties and equation of state of nuclear matter are assesed and analyzed in the framework 
of the proposed effective interaction model. Special attention is dedicated to the study of 
>> ■ the contribution of the components of /3-stable nuclear matter to the entropy per particle, a 
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1 Introduction 

The equation of state (EOS) of hot nuclear matter determines the structure inside a supernova p] 
and a hot neutron star [21 El S] and affects the state of matter, such as its chemical composition. 
In addition, the equation of state plays important roles for the study of the supernova explosion, 
as well as on determining the evolution of a neutron star at the birth stage. The profiles of a 
neutron star as the density, temperature and proton fraction during the cooling, which affect the 
reaction rate of neutrino process inside the star, are determined through the equation of state. 

There exist many calculations for hot nuclear matter with applications to the properties of hot 
neutron stars and supernova [HElIallllElIHlITllHlEl HOl HH [121 H31 HH HSl [161 HZl HH [19l [20l [2H 
12211231211231201123123123^ G. Baym et al. provided 

an EOS of neutron matter [5] and Bethe et al. [B] an EOS for the gravitational collapse of stars. 
Friedman and Pandharipande [7] performed variational calculations of the equation of state of 
hot and cold, nuclear and neutron matter. Lattimer and Swesty carried out calculations of the 
EOS for stellar collapse, using the compressible liquid-drop model for nuclei jBj. M. Prakash et 
al. [2] investigated the structure of neutron stars shortly after their birth, by applying various 
nuclear models. Takatsuka et al. [3 \TU\ have performed detailed calculations for supernova matter, 
within the framework of finite temperature Hartree-Fock approach, with effective nucleon-nucleon 
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interaction. Recently Das et al. [TT] have calculated the EOS of dense supernova matter within 
the finite temperature Brueckner Goldstone approach with effective two-body Sussex interaction. 

The present work is based on the previous work of Prakash et al. [2J. More specifically, in 
order to study the properties and the EOS of hot nuclear matter, we apply a momentum-dependent 
effective interaction model (MDIM), which is able to reproduce the results of more microscopic 
calculations of dense matter at zero temperature and which can be extended to finite temperature 

[a sum]. 

The aim of this work is to apply a momentum-dependent interaction model for the study of 
hot nuclear matter EOS under /5-equilibrium. The present model has the additional property, 
compared to the previous ones, that the temperature affects not only the kinetic part of the en- 
ergy density, but also influences the interaction part of the energy density as well. In that way, 
we are able to study simultaneously thermal effects not only on the kinetic part of the symmetry 
energy and symmetry free energy, but also on the interaction part of the above quantities. This 
is important in the sense that the density dependent behavior of the symmetry energy and sym- 
metry free energy influence strongly the values of the proton fraction and as a consequence the 
composition of hot /3-stable nuclear matter, under consideration. 

Using the above method , we will show that the thermal energy (and also the related quantities) 
depend sensitively on the momentum dependence of the nuclear interaction. We concentrate our 
study on the properties of hot nuclear matter in the density range < n < Quq (where uq = 0.16 
fm -3 is the saturation density) and temperature range < T < 30 MeV, taking into account that 
nuclear matter consists of neutrons, protons, electrons and muons with their relative concentrations 
determined from the conditions of charge neutrality and equilibrium under /3-decay process in the 
absence of neutrino trapping. 

The article is organized as follows. In Sec. II the model and relative formulas are discussed 
and analyzed. Results are reported and discussed in Sec. Ill, whereas the summary of the work 
is given in Sec. IV. 

2 The model 

We start by outlining the momentum dependent interaction model, then we define the thermo- 
dynamic quantities of nuclear matter and finally we analyze the /3-equilibrium conditions, the 
contribution on pressure and energy of leptons and the total equation of state of nuclear matter. 

2.1 Momentum dependent interaction model 

The schematic potential model, employed here, is designed to reproduce the results of the mi- 
croscopic calculations of both nuclear and neutron-rich matter at zero temperature and can be 
extended to finite temperature [2J. The energy density of the asymmetric nuclear matter (ANM) 
is given by the relation 

e(n n , n p , T) = e kin (n n , T) + e p kin (n p , T) + V int (n n} n p , T), (1) 

where n n (n p ) is the neutron (proton) density and the total baryon density is n = n n + n p . The 
contributions of the kinetic parts are 

r d 3 k Ti 2 k 2 

e n kin (n n , T) + e p kin (n p , T) = 2 J j^y 3 ^ (/«K, k, T) + f p (n p , k, T)) , (2) 
where f T , (for r = n,p) is the Fermi-Dirac distribution function with the form 

f T (n T ,k,T) = 



1 + exp 



k,T) - fi T (n T ,T) 



(3) 
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The nucleon density n T is evaluated from the following integral 

'e T (n T ,k,T) - / u T (n r ,T)~ 



d 3 k r d 3 k 

f T {n TJ k,T) = 2 



(2tt 



(27T) 



1 + exp 



T 



-i 



(4) 



In Eq. d3J), e T (n T , k,T) is the single particle energy (SPE) and // T (n r ,T) stands for the chemical 
potential of each species. The SPE has the form 



e T (n T , k,T) 



2 7.2 



K'k 

2m 



+ £y r (n r , fc,T), 



(5) 



where the single particle potential U T (n T , k,T) is obtained by the functional derivative of the 
interaction part of the energy density with respect to the distribution function f T . Including the 
effect of finite-range forces between nucleons, to avoid acausal behavior at high densities, the 
potential contribution is parameterized as follows [2] 



V int {n n) n P) T) 



-Ann 



x )P 



u 



\Bn 


|-(i + x 3 )/ 2 


u a+l 


1 + lB' 




u a - 1 



+ u £ 



where 



i=l,2 



Jr 



C% ( J7n + J7p) H - I (J7n ~~ <Jp) 



d 3 k 



g(k,Ai)f T (n T ,k,T). 



(6) 



(7) 



(2vr) 3 - 

In Eq. (EJ), I is the asymmetry parameter (/ = (n n — n p )/n) and u = u/uq, with no denoting 
the equilibrium symmetric nuclear matter density, hq = 0.16 fm~ 3 . The asymmetry parameter / 
is related to the proton fraction Y p by the equation I = (1 — 2Y P ). The parameters A, B, a, C\, 
C*2 and B', which appear in the description of symmetric nuclear matter, are determined in order 
that E(n = n ) — mc 2 = —16 MeV, n = 0.16 fm~ 3 , and the incompressibility to be K = 240 
MeV. The additional parameters Xq, £3, Z\, and Z 2 , which are used to determine the properties 
of asymmetric nuclear matter, are treated as parameters constrained by empirical knowledge [2]. 
The parameterizations used in the present model have only a modest microscopic foundation. 
Nonetheless, they have the merit of being able to closely approximate more physically motivated 
calculations as presented in Fig. 1. More precisely, in Fig. 1 we compare the energy per baryon 
(for symmetric nuclear matter (Fig. la) and pure neutron matter (Fig. lb)) calculated by the 
present schematic model (MDIM), with those of existent, state of the art calculations by Wiringa 
et al. [JT] and Pandharipande et al. [4~2] . 

The first two terms of the right-hand side of Eq. arise from local contact nuclear interaction 
which lead to power density contributions as in the standard Skyrme equation of state. The 
first one concerns attractive interaction while the second one is repulsive, and both are assumed 
to be temperature independent. The third term describes the effects of finite range interactions 
according to the chosen function g(k, Aj), and is the temperature dependent part of the interaction. 
This interaction is attractive and important at low momentum, but it weakens and disappears at 
very high momentum. The function g(k,Ai), suitably chosen to simulate finite range effects, has 
the following form [2] 

'Jl 



g{k,A) 



1 



-1 



where the finite range parameters are Ai 
at the saturation point no- 



l.bkp and A2 = 3k F and k° F is the Fermi momentum 
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The main origin of the momentum dependence in Brueckner theory is the nonlocality of the 
exchange interaction. Following the discussion of Bertsch et al. [13] a single-particle potential 
U (n) which depends only on the baryon density is oversimplified. What is more, it is well known 
that nuclear interaction has strong exchange effects which give rise to a momentum dependence 
in the single-particle potential and as a consequence has an effect on the energy density func- 
tional. The question here is how best to parameterize the momentum dependence in modelling 
the potential U(n, k). A promising approach might be to adopt the relativistic mean field model, 



where U(n, k) = U v n + 



UsTl 



The above potential exhibits a strong momentum dependence 



for small k which diminishes to zero at high momentum. In order to perform extensive studies in 
heavy ion collision studies, Gale et al. [H] have proposed the following parametrization for the 
momentum part of the single-particle 



U(n,k) ~ C 



n 



n l + (k- (k'» 2 /A 2 ' 



This has a proper fall-off at high k and Galilean invariance is assured by measuring k with respect 
to the average of the particles in the neighborhood, (k'). For static nuclear matter we have 
(k') = 0. 

The present model, which is a generalization of that proposed by Gale et al. [B], has been 
successfully applied in heavy ion collisions and astrophysical studies over the years El HI E3 EEl 

[331 mm]. 

In order to clarify the relative contribution of the three terms of the potential energy density 
mentioned above, we plot them as a function of the baryon density in Fig. 2a. In this figure we 
have that 



Va 
V b 

V c 



-An n 



Bn 



f-<5 + *> 



u 



u 



cr+l 



l + \B> §-(§ + x 3 )P 



u 



a-l 



(9) 



8=1,2 



c 1 [j: 1 + j;) + 



{Ci - 8Zi) 



As mentioned above, the first term Va corresponds to an attractive interaction, whereas the second 
one Vb corresponds to repulsive interaction and dominates for high values of n (n > 0.6 fm~ 3 ). 
Both of these terms are temperature independent. The third term Vc (is plotted for T = 0) 
contains the momentum dependent part of the interaction, corresponds to attractive interaction, 
and its main contribution is to compete with the repulsive interaction of Vb for high values of n 
and as a consequence avoid acausal behavior of the EOS at high densities. The term Vc consists of 
two finite range terms, one corresponding to a long-range attraction and the other to a short-range 
repulsion. 

Thermal effects on the momentum dependent term Vc are displayed in Fig. 2b. The contri- 
bution of Vc is plotted for various values of T. It is therefore concluded that thermal effects are 
more pronounced for high values of T (T > 10 MeV), leading to a less attractive contribution. 
More precisely, we find that for small values of n (i.e. n = 0.15 fm~ 3 ) Vc increases (compared to 
the cold case T = 0) 3% — 20% for T = 10 — 30. For higher values of n the increase is even less. 

An additional test for the present model is to compare the single particle potential U T (n T , k, T) 
(or U T (n, I, k, T)) originated from the present version of the momentum dependent interaction with 
other calculations. The single particle potential U T (n, I, k, T) (protons or neutrons), obtained from 
the functional derivative of the interaction part of the energy density (Eq. <^j) with respect to the 
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distribution function f T , has the general form [10] 

U T (n, I, k, T) = U?(n, I) + U?(n, I) + Uf D (n, I, k, T). (10) 

It is of interest to see that the single particle potentials are separated into two parts. The first 
one, U^(n,I) + U^(n,I) is a function only of the baryon density n and the isospin asymmetry 
parameter I. The second one, U^ D (n, I, k, T) has an additional dependence on T and k. Actually, 
U^ ID {n, I, k, T) is mainly responsible for the trend of the effective mass and also the effective mass 
splitting. Additionally, it is connected with the effect of the temperature on the interacting part 
of the energy density [10] ■ 

The single-particle potential in symmetric nuclear matter has been calculated microscopically 
for several Hamiltonians by Wiringa [47]. These Hamiltonians include nucleon-nucleon potentials 
fit to scattering data and three nucleon potentials fit to binding energies of few-body nuclei and 
saturation properties of nuclear matter. The potential was parameterized using the ansatz 

U(n,k) = a(n)+ (11) 

where the density dependent parameters a(n), (3{n) and A(n), for three types Hamiltonians, are 
listed in Table I of Ref. [47] . 

Furthermore, the single-particle potential has been derived by Li et al. [48]. The derivation is 
based upon the Bonn meson-exchange model for the nucleon-nucleon interaction and the Dirac- 
Brueckner approach for nuclear matter. The potential, named DBHF, has been parameterized as 
following 

U(n, k) = an + /?n 7 + 5 \n 2 (e(hck) 2 + l)rf . (12) 

The parameters a, /3, 7, 5, e and a are listed in Table I of Ref. [38]. 

A comparison with the results of UV14+TNI, UV14+UVII, AV14+UVII and DBHF interac- 
tions show that (see Fig. 3) U(n, k) (for T = 0) obtained from the present model is very reasonable 
for at least up to the value k = 3 fm _1 . The agreement is not so good for high values of k, es- 
pecially compared with the DBHF interaction, but as has been pointed out by Li et al. [H], the 
parametrization in Eq. [12] is bad for large k since it continues to grow with increasing k, while the 
exact potential becomes independent of k for large momenta. In conclusion, the present results 
show that the momentum dependent interaction model, which has been applied in the present 
work, provides a reliable representation of U(n, k) for a wide range of density and momentum. 

The energy density of asymmetric nuclear matter at density n and temperature T, in a good 
approximation, is expressed as 

e(n, T, /) = e(n, T,I = 0) + e sym (n, T, 7), (13) 

where 

e sym (n, T, I) = nI 2 E™ m {n, T) = nl 2 (J5j*>, T) + E^ m (n, T)) . (14) 

In Eq. (fl4"j) the nuclear symmetry energy E%jf m (n,T) is separated in two parts corresponding to 
the kinetic contribution E^™(n,T) and the interaction one E'%£ n {n,T). 

From Eqs. ffTBT) and (|14|) and setting / = 1, we find that the nuclear symmetry energy 

Elfjn, T) = i (e(n, T, I = 1) - e(n, T, 1 = 0)). (15) 

Thus, from Eq. ( ITBT) and by a suitable choice of the parameters Xq, £3, Z\ and Z 2 , we can obtain 
various forms for the density dependence of the symmetry energy El?f m (n,T). 
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It is well known that the need to explore different forms for El°f m (n,T) stems from the uncer- 
tain behavior at high density [2j. The high-density behavior of symmetry energy is the least known 
property of dense matter [HH [50j [51], with different nuclear models giving contradictory predic- 
tions. Thus, in relativistic mean field (RMF) models, the symmetry energy increases strongly 
with the density of nuclear matter [52] , while in many realistic potential models of nuclear matter 
in the variational approach [7J [53], the symmetry energy saturates and then bends over at higher 
densities. 

Recently, the density dependence of the symmetry energy in the equation of state of isospin 
asymmetric nuclear matter has been studied using isoscaling of the fragment yields and the an- 
tisymmetrized molecular dynamic calculation [54|. It was observed that the experimental data 
at low densities are consistent with the form of symmetry energy, E sym (u) ~ 31.6m 0,69 , in close 
agreement with those predicted by the results of variational many-body calculations. In Ref. [51] 
it was suggested also that the heavy ion studies favor a dependence of the form E sym {u) ~ 31.6m 7 , 
where 7 = 0.6 — 1.05. This constrains the form of the density dependence of the symmetry energy 
at higher densities, ruling out an extremely "stiff" and "soft" dependence [51]. 

Additionally, Chen et al. [55] also showed, using the isospin dependent Boltzmann-Uehling- 
Uhlenbeck transport model calculations, that a stiff density dependence of the symmetry energy 
parameterized as, E sym (u) ~ 31. 6u 105 explains well the isospin diffusion data [56J from NSCL- 
MSU (National Superconducting Cyclotron Laboratory at Michigan State University). 

In this paper, since we are interested mainly in the study of thermal effects on the nuclear 
symmetry energy, we choose a specific form for it, enabling us to reproduce accurately the results 
of many other theoretical studies [571 EH]- In Ref. [57] the authors carried out a systematic analysis 
of the nuclear symmetry energy in the formalism of the relativistic Dirac-Brueckner-Hartree-Fock 
approach, using the Bonn one-boson-exchange potential. In a very recent work |58j, the authors 
applied a similar method as in Ref. [57] for the microscopic predictions of the equation of state in 
asymmetric nuclear matter. In that case E sym (u) is obtained employing the simple parametrization 
E sym (u) = Cu 1 with 7 = 0.7— 1.0 and C ~ 32 MeV. The authors conclude that a value of 7 close 
to 0.8 gives a reasonable description of their predictions, although the use of different functions in 
different density regions may be best for an optimal fit [58] ■ The results of Refs. [57] 158] are well 
reproduced by parameterizing the nuclear symmetry energy according to the formula 

EllUn, T = 0) = 13v^ + lTF^ . (16) 

Kinetic Interaction 

For the function F(u), which parametrizes the interaction part of the symmetry energy, we apply 
the following form 

F(u)=u. (17) 

The parameters x , x 3 , Z 1 and Z 2 are chosen so that Eq. (TT5]) . for T = 0, reproduces the results 
of Eq. (TTBl) for the function F{u) = u. 

In one of our previous paper [5H] , the potential part of the symmetry energy has been parame- 
terized in the generalized form, F{u) = u c , and the obtained nuclear equations of state are applied 
to the systematic study of the global properties of a neutron star (masses, radii and composition). 
We obtained a linear relation between the parameter c and the radius and the maximum mass of 
the neutron star [59J. Additionally, we found that a linear relation between the radius and the 
derivative of the symmetry energy near the saturation density n holds [59J . 

It is worthwhile to point out that the above parametrization of the interacting part of the 
nuclear symmetry energy is used extensively for the study of neutron star properties [21 [60] , as well 
as for the study of the collisions of neutron- rich heavy ions at intermediate energies [611 [62] • For a 
very recent review of the applications of the proposed momentum dependent effective interaction 
model and its specific parameterizations, see Ref. [1] (and references therein). 
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2.2 Thermodynamic description of hot nuclear matter 

In order to study the properties of nuclear matter at finite temperature, we need to introduce 
the Helmholtz free energy F. The differential of the total free energy F to t (the total free energy 
of baryons contained in volume V) and total internal energy E tot (the total internal energy of 
baryons contained in volume V) are given by 



dF tot = -S tot dT -PdV + J2 (18) 

i 

dE tot = TdS tot -PdV + J2 fMdNt (19) 

i 

where S to t is the total entropy of the baryons, while //j and iV, are the chemical potential and the 
number of particles of each species respectively. 

It is easy to prove that the free energy per particle F is written as 



F(n, T, I) = E(n, T, I) - TS{n, I, T). (20) 

In Eq. (1201) . E is the internal energy per particle, E = e/n, and S is the entropy per particle, 
S = s/n. From Eq. (T2"0l is also concluded that for T = 0, the free energy F and the internal 
energy E coincide. 

The entropy density s has the same functional form as that of a non interacting gas system, 
given by the equation 

/d h 
j— 3 [f T In f T + (1 - f T ) ln(l - f T )] . (21) 

The total internal energy E tot is useful for studying isentropic processes. In that description of a 
thermodynamic system, the pressure and the chemical potential are defined as follows 



fdE tot \ Jde/n\ ( dE tot \ ( de\ 

P = ~[~dV r ) =U ' ^ = I ~dN~ ) = [d^- • (22) 

Now we are going to study the properties and the EOS of nuclear matter by considering 
an isothermal process. In that, the pressure and the chemical potential are connected with the 
derivative of the total free energy F tot . More specifically, they are defined as follows 

p = _ ( dF **\ = n 2 ( d f/ n ) . = ( 9F toA = (M\ 

where / is the free energy density. 

The pressure P can be calculated also from equations 



(23) 



VP = TS tat -E tot + ^2mN i , or P = Ts-e + J2vini. (24) 

i i 

It is also possible to calculate the entropy per particle S(n, T) by differentiating the free energy 
density / with respect to the temperature 

S(n,T) = -( d -m . (25) 

V U1 / V,Ni 

The comparison of the two entropies, that is from Eqs. (1211 and ( 1251) . provides a test of the 
approximation used in the present work. 
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It is easy to show by applying Eq. fl2"31 that (see for a proof [HO] as well as 



du ) Y _ T \9Y f 



P/ n.T 

dF\ 



\ UI P/ n,T 

(OF \ 
qyJ ■ ( 26 ) 

We can define the symmetry free energy per particle F sym (n,T) by the following parabolic ap- 
proximation (see also [32| [33]) 

F(n, T, I) = F(n, T, / = 0) + I 2 F sym (n, T) = F(n, T, I = 0) + (1 - 2Y p ) 2 F sym (n, T), (27) 

where 

F sym (n, T) = F(n, T, I = 1) - F(n, T, I = 0). (28) 

It is worthwhile to notice that the above approximation is not valid from the beginning, but one 
needs to check the validity of the parabolic law in the present model before using it. As we see 
later, that law is well satisfied as well as the parabolic law holding for the energy. 
Now, by applying Eq. (I2TI) in Eq. (121))) . we obtain the key relation 

A = Mn - M P = 4(1 - 2Y p )F sym {n, T). (29) 
The above equation is similar to that obtained for cold nuclear matter by replacing E sym (n) with 



2.3 /^-equilibrium, leptons contribution and equation of state 

Stable high density nuclear matter must be in chemical equilibrium for all types of reactions, 
including the weak interactions, while f3 decay and electron capture take place simultaneously 

n — > p + + u e , p + e~ — ► n + v e . (30) 

Both types of reactions change the electron per nucleon fraction, Y e and thus affect the equation 
of state. Here, we assume that neutrinos generated in those reactions have left the system. The 
absence of neutrino-trapping has a dramatic effect on the equation of state and mainly induces a 
significant change on the values of the proton fraction Y p [91 [10] . The absence of neutrinos implies 
that 

p, = H n - ll p = fi e . (31) 

When the energy of electrons is large enough (i.e. greater than the muon mass), it is energet- 
ically favorable for the electrons to convert to muons 

e~ — > n~ + Vy, + v e . (32) 

Denoting the muon chemical potential by the chemical equilibrium established by the above 
process and its inverse is given by 

A* e = /V 

Taking into account that the threshold for muons occurs for /z M = m^c 2 ~ 105.7 MeV, one may 
expect muons to appear roughly at nuclear density n = 0.16 fm -3 . 
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Thus, in total, we consider that nuclear matter contains neutrons, protons, electrons, and 
muons. They are in a /3-equilibrium, where the following relations hold 



Furthermore, they obey the charge neutrality condition i.e. 

n p = n e + rifj,. 

The leptons (electrons and muons) density is given by the expression 



n, 



dk 



1 + exp 



^Jh 2 k 2 c 2 +m 2 c A — f±i 
T 



(33) 



(34) 



(35) 



One can solve self-consistently Eqs. (j29]) . (l33]) . (l33j) and (I3"5"j) in order to calculate the proton fraction 
Y p , the lepton fractions Y e and Y^, as well as the electron chemical potential \x e as a function of 
the baryon density n, for various values of the temperature T. 

The next step is to calculate the energy and pressure of leptons given by the following formulae 



ei{n h T) 



(27T) 



h 2 k 2 c 2 + m 2 c A dk 



1 + exp 



T 



k 2 dk 



yJh 2 k 2 c 2 +m 2 c A — [it 



(36) 



(37) 



3 (27T) 3 J ^ k 2 c 2 + m 2 c 4 1 + 

The chemical potentials of electrons and muons are equal and according to Eqs. (1291 and (1331) are 

fJ-e = AV = Hp ~ = 4(1 - 2Y p (n,T)) F sym (n,T) = AI(n,T)F sym (n,T). (38) 

The equation of state of hot nuclear matter in /3-equilibrium (considering that it consists of 
neutrons, protons, electrons and muons) can be obtained by calculating the total energy density 
e to t as well as the total pressure P tat . The total energy density is given by 



e tot (n, T, I) = e b (n, T, I) + e K n > T > J )> 



(39) 



-e,fj, 



where e b (n,T,I) and ei(n,T,I) are the contributions of baryons and leptons respectively. The 
total pressure is 

P tot (n,TJ)=P b (n,T,I) + £ Pi(n,T,I) 7 (40) 

l=e,fi 



where P b (n,T,I) is the contribution of the baryons (see Eq. (1241) ) i.e. 

P b (n,T,I) =T s T (n,T,I)+ n T /i T (n,T, I) - e b (n,T, I), 



(41) 



r=p,n 



T=n,p 



while Pi(n, T, I) is the contribution of the leptons (see Eq. (1371) ). From Eqs. (1391 and (1401 we can 
construct the isothermal curves for energy and pressure and finally derive the isothermal behavior 
of the equation of state of hot nuclear matter under /5-equilibrium. 
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3 Results and Discussion 



The schematic potential model, which has been applied in the present work, has been designed 
to reproduce the results of the more microscopic calculations of both nuclear and neutron-rich 
matter up to high values of baryon density (see Fig. 1). The behavior of the high density EOS 
is of great significance to the determination of hot protoneutron stars and cold neutron stars 
structure. The model has the additional advantage that with the appropriate parametrization, is 
able to reproduce different forms of the density dependence of the nuclear symmetry energy. 

In view of the above discussion we calculate the equation of state of hot asymmetric nuclear 
matter by applying a momentum dependent effective interaction model describing the baryons 
interaction. We consider that nuclear matter contains neutrons, protons, electrons and muons 
under /3-equilibrium and charge neutrality. The key quantities in our calculations are the proton 
fraction Y p and also the asymmetry free energy defined in Eq. (1281) . 

In order to check the validity of the parabolic approximation (|2"Tj) . we plot in Fig. 4 the 
difference F(n,T,I = 1) — F(n,T,I = 0) as a function of (1 — 2Y P ) 2 at temperature T = 10 
and T = 30 MeV for three baryon densities, i.e., n = 0.2 fm -3 , n = 0.3 fm -3 , and n = 0.4 
fm -3 . It is seen that in a good approximation an almost linear relation holds between F(n, T,I — 
1) — F(n,T,I = 0) and (1 — 2Y P ) 2 . A similar behavior of F sym (n,T) is found by Xu et al. [33] . 
applying an isospin and momentum dependent interaction model. 

It is worth to present the calculation recipe of our work. The outline of our approach is the 
following: For a fixed baryon density n, temperature T, and asymmetry parameter J, Eq. (J4J) may 
be solved iteratively in order to calculate the quantity 

»(n r ,D = Mn " r) ; R( "" r) , (42) 

where 

U T {n T , T) = U T (n T , k, T) - U T (n T , k). (43) 

Knowledge of r] T (n, T) allows the evaluation of U T (n T , T), which then may be employed to infer 
the chemical potential from 

Li T (n T ,T) = T VT (n T ,T) + U(n T ,T), (44) 

required as an input for the calculation of the Fermi-Dirac distribution function f T (n T , k,T). The 
knowledge of f T (n T ,k,T) permits the calculation of the bulk quantities of asymmetric nuclear 
matter. 

F sy m(n, T), for various values of the temperature T, was derived with a least-squares fit to the 
numerical values according to Eq. (|28|) and has the form 

F sym {u; T = 0) = 13w 2/3 + 17m 

F sym (u;T = 5) = 3.653 + 28.018m -1.512u 2 + 0.185m 3 -O.OOlu 4 , 

F sym (u; T = 10) = 5.995 + 26.157m - 0.827u 2 + 0.068u 3 - 0.002u 4 , 

F sy m(u; T = 20) = 13.200 + 21.267m + 0.800m 2 - 0.193m 3 + 0.014m 4 , 

F sym (u;T = 30) = 21.087 + 17.626u + 1.645m 2 - 0.289m 3 + 0.018m 4 . (45) 

where the case with T = 0, is included as well. In that case F sym coincides with E sym . 

In Fig. 5 we present the behavior of the free energy, corresponding to hot /3-stable nuclear 
matter, as a function of the baryon density n, for various values of the temperature T. It is 
obvious that the thermal effects are more pronounced for low values of the density n. 
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In Fig. 6 we plot the calculated free energy for symmetric nuclear matter and pure neutron 
matter of the proposed momentum dependent interaction model in comparison with the values 
of the free energy calculated by Friedman and Pandharipande (FP model) [7] . In the FP model 
the equation of state of hot and cold nuclear and neutron matter has been calculated in the 
framework of a variational calculation, where a realistic nuclear interaction containing two- and 
three-nucleon body nucleon-nucleon interaction has been used. In the case of symmetric nuclear 
matter the results of the two models are very similar up to values n = 0.4 — 0.5 fm~ 3 depending 
on the values of T. The above agreement is expected, in the sense that a part of the parameters of 
the model applied in the present work are determined from constraints provided by the empirical 
properties of symmetric nuclear matter at the equilibrium density hq = 0.16 fm~ 3 . 

However, there is an obvious disagreement in the case of pure nuclear matter, where in the 
two models the free energy exhibits a different trend, especially for higher values of n. The above 
disagreement will be explained below. 

In Fig. 7 we display the internal energy per particle E(n,T) = e(n,T)/n given by Eq. ([T]) for 
various values of temperature. Thermal effects, as expected, are more pronounced for low values 
of the baryon density n and less important for high values of n. 

In Fig. 8 we display the internal energy of symmetric nuclear matter and pure neutron matter, 
for T = 0, calculated by the MDIM and FP models. In accordance with the case of the free 
energy, there is a very good agreement in symmetric nuclear matter, but an obvious disagreement 
is exhibited in pure neutron matter. The explanation of the agreement in the first case is the same 
as in the case of the free energy. The disagreement is due to the completely different behavior of 
the two models of the nuclear symmetry energy, presented in Fig. 9. In our model the parameters 
xo, xs, Z\ and Z 2 chosen so that Eq. (ITS]) , for T = 0, reproduce the results of Eq. ([16]) for the 
function F[u) = u. Consequently, E sym (n) shows an increasing trend shown in Fig. 9. In contrast, 
in the FP model, E sym (n) is a slightly increasing function of n for low n and then a decreasing 
function of n for n > 0.5 fm™ 3 . 

In addition, we plot the nuclear symmetry energy extracted from experimental results and pre- 
sented in Ref. [51], where E sym {u) is parameterized according to the relation E sym {u) « 31.6m 0-69 
as well as experimental results extracted from Ref. [55] , where E sym {u) is given by E sym {u) ~ 
31.6m 105 . The important point to be noted is that both cases clearly favor a stiff density de- 
pendence of the symmetry energy at higher densities, ruling out the very stiff and very soft 
predictions. These results can thus be employed to constrain the form of the density dependence 
of the symmetry energy at supranormal densities relevant for neutron star studies [SI]. In the 
same figure the theoretical predictions of Ref. |5Ej are presented, where E sym {u) is parameterized 
by E sym (u) « 32m 08 . 

The results of Ref. [SI] are in a good agrement with the present model up to n = 0.3 fm -3 
while the theoretical predictions of Ref. [58] are very close to the present model up to very high 
values of the baryon density n. 

However, our motivation, here, is not to perform a systematic comparison of various models, 
but we would like just to present the similarities and the deviations existing between the models. 
The deviations, concerning the symmetry energy behavior of the two models (MDIM and FP 
model) are well reflected on the behavior of the free energy and internal energy of pure neutron 
matter as shown in Figs. 6 and 8. 

In Fig. 10 we plot the thermal energy per particle 

E t hermal(n, T) = E(n, T) - E(n, T = 0), 

of /3-stable matter as a function of the baryon density n for various values of temperature T. The 
most striking feature of i?fft erma /(n, T) is that for small values of T, the thermal contribution to 
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the internal energy is almost independent of the density n. For high values of T the situation is 
different and E t h e rmai(n,T), for fixed values of T, is a decreasing function of the density n. 

Ethermai(n, T) can be decomposed to separate contributions of the kinetic and potential energies 
as follows: 

Ethermaiin^T) = E^ rmal (n,T) + Ef hermal (n, T) . 

We find that for small values of the baryon density (i.e. n = 0.2 fm -3 ) the contribution, to 
E t hermai{n,T) of the potential energy 

E ?hermai( n , T ) is 20% - 10% for T = 5 - 30 MeV. For 
medium values of n (i.e. n = 0.4 fm -3 ) is 43% — 20% for T = 5 — 30 MeV and for higher values 
(i.e. n = 0.6 fm -3 ) is 70% — 30% for T = 5 — 30 MeV. Hence, it is concluded that the potential part 
of the energy (as a result of the momentum dependence of the interaction) contributes significantly 
to the thermal energy, mainly for small values of T (for fixed values of n) and for large values of 
n (for fixed values of T) . 

At this point, it is worthwhile to compare the results for the pressure obtained by applying 
Eqs. ( |24"|) and (|23|) . Thus, in Fig. 11 we plot P of asymmetric nuclear matter for Y p = 0.1 and 0.3 
at T = 10 and 30 MeV. The full lines give the results calculated from Eq. (T24|) . while the squares 
give results obtained by differentiating F(n,T) (Eq. (!23l ). The two calculations for the pressure 
are in excellent agreement. This agreement provides a test of the calculations performed in the 
present model. 

It is of interest also to study the effect of the temperature on the baryon pressure defined by 
equation (1241) . A related quantity is the thermal pressure Pthermai(n,T) defined as: 

Pthermal(n, T) = P(n, T) - P(n, T = 0). 

T) as a function of n, for various values of T is seen in Fig. 12. P f / ierma /(n, T), in all of 
the cases, is an increasing function of the baryon density. 

The proton fraction affects the reaction rate of neutrino process inside that star. If a neutron 
star has a large proton fraction, the cooling rate may drastically change through the high neutrino 
emissivity due to the direct Urea process. This process can occur if the proton fraction in the 
matter of a cold neutron star exceeds the critical value of 0.11-0.15 and would lead to the rapid 
cooling of the neutron star. Thus, it is important to calculate the proton fraction as a function of 
the baryon density and investigate the temperature effects on that. 

Fig. 13 displays the fractions of protons, electrons and muons as functions of the density, for 
various values of T. The proton fraction is an increasing function of T and this effect is more 
pronounced for T > 10 MeV. The proton fraction Y p was derived also with a least-squares fit to the 
numerical results obtained from our calculations, leading to the following relations (for n > 0.15 
fm -3 ). 



Y p {n-T~- 


= 0) 


= -0.050 + 0.633n-0.521n 2 + 0.184n 3 , 


Y p (n;T-- 


= 5) 


= -0.046 + 0.625n-0.514n 2 + 0.179n 3 , 


Y p (n;T = 


10) 


= -0.032 + 0.570n- 0.436n 2 + 0.139n 3 , 


Y p (n;T = 


20) 


= 0.021 + 0.378n-0.163n 2 + 0.004n 3 , 


Y p (n;T = 


30) 


= 0.109 - 0.062n + 0.908n 2 - 1.270n 3 + 0.580n 4 



In Fig. 14 we plot the Fermi distribution function / p n (n, T) both for neutrons and protons for 
various values of T. We observe that the diffuseness of f p {n,T) is larger than that of f n (n,T). 
We give an explanation (see also [9]): the ratio of T to the Fermi kinetic energy epi is a measure 
of the thermal effect. Thus by comparing the two ratios we have (see also Appendix) 

( ( T /^p) ^ _ h n2/ k p2 _ y3 /2/y3/2 
\ frr I \ I F "'P n / 1 p ■ 
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But, due to 7 n > Y p , we conclude that we expect the diffuseness to be larger for the proton 
distribution than for the neutron one, depending of course on the specific value of the ratio Y p /Y n . 
As we will see later, this fact influences the values of the contributions of protons and neutrons to 
the total entropy per particle. The entropy, according to relation (12TI) . is an increasing function of 
the diffuseness. Thus, the contribution of each species on the total value of the entropy depends 
strongly on the diffuseness of the corresponding Fermi distribution function. 

We can provide a second test, concerning the accuracy of the present calculations, by comparing 
the entropy per baryon calculated by applying Eqs. (I2T]) and (|25|) . Thus, in Fig. 15 we plot S of 
asymmetric nuclear matter with Y p = 0.2 at T = 10, 20, 30 MeV. The full lines give the entropy 
calculated from Eq. fl2Tj) . while the squares give results obtained by differentiating F(n,T) with 
respect to T (Eq. §25§). It is obvious that there is again a very good agreement of the results, 
especially for low values of T and high values of n. 

In Fig. 16 we plot the contribution of the proton S p , the neutron S n and the total entropy 
per baryon S. It is obvious that there is a strong effect of T on the values of the entropies 
mainly for low values of the density. The main part of the contribution comes from neutrons, 
whereas the contribution of protons is three times less. It is worthwhile to notice that, in spite 
of Y p ~ (1/20 — 3/10)F n , the approximate relation S p ~ (1/4 — 3/7)S n holds. This feature is 
understood by the previous discussion that f p (n, T) is diffused more broadly than f n (n, T), so the 
larger the diffuseness, the larger is the entropy contribution (see also [H]). 

In Fig. 17 we plot the contribution of the electronic S e , the muonic S M and the total (leptonic) 
Si to the entropy per baryon. The contribution to the entropy, of S e depends slightly on the 
density, for fixed values of T. Our present results are very close to those found by Onsi et al. [19] , 
where they employed the analytical approximate formula for the electron entropy density s e 

s e = lj^- 3 T ' A*. = nc{^ 2 Y e nf/\ (47) 

According to the above formula, the contribution of electrons to the entropy per baryon has the 
form 

/y2\l/3 

S e = a e /n~\-±.\ T. (48) 

/ y2 \ 1/3 

The quantity f- 2 -) is a function slightly dependent on the density n, so that for a fixed value 
of T the contribution S e is almost constant. The muonic contribution to the entropy, for fixed T, 
increases slightly SIS db function of the density. 

In Fig. 18 we present the EOS of the /3-stable hot nuclear matter by taking into account and 
analyzing the contribution to the total pressure of each component. The main contribution to the 
total pressure originates from the baryons, while the contribution of the leptons is about a few 
percent compared to P^. It is worthwhile to notice that thermal effects are not important for the 
calculation of P e , but only for P M , especially for small values of n (n < 0.4 fm~ 3 ). We found that 
thermal effects produce a slightly stiffer equation of state with respect to the case of cold nuclear 
matter. The above EOS can be applied to the evaluation of the bulk properties of hot neutron 
stars (mass and radius). 

The study of hot nuclear matter in the absence of neutrino trapping is the first step to study 
the properties of hot neutron stars and supernova matter. Next, one can study the more realistic 
case of neutrino-trapped matter in /3-equilibrium. In this case, the /3-equilibrium conditions in 
matter are altered from the case in which neutrinos have left the system and thus the composition 
of matter is affected. The proton fraction increases dramatically and influences significantly the 
properties of nuclear matter. Such a work is in progress. 
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4 Summary 



The evaluation of the equation of state of hot nuclear matter is an important problem in nuclear 
physics and astrophysics. EOS is the basis ingredient for the study of the supernova explosion as 
well as on determining the properties of hot neutron stars. The motivation of the present work is 
to apply a momentum-dependent interaction model for the study of the hot nuclear matter EOS 
under /3-equilibrium in order to be able to study simultaneously thermal effects, not only on the 
kinetic part of the symmetry energy and symmetry free energy, but also on the interaction part 
of the above quantities as well. We calculate the proton fraction, as well as the lepton fractions, 
by applying the constraints for chemical equilibrium and charge neutrality. The free energy, the 
internal energy and also the pressure are calculated as functions of baryon density and for various 
values of temperature. We also concentrate on the evaluation of thermal effects on the internal 
energy and baryon pressure. Special attention is dedicated to the study of the contribution of the 
components of /3-stable nuclear matter on the entropy per particle, a quantity of great interest in 
the study of structure and collapse of supernova. We present and analyze the contribution of each 
component. Finally, we present the EOS of /3-stable hot nuclear matter, by taking into account 
and analyzing the contributions to the total pressure of each component. The above EOS can be 
applied to the evaluation of the gross properties of hot neutron stars i.e. mass and radius, (work 
in progress). 
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Appendix 

The energy density of baryons (Eq. (pQ)), at T = 0, is given by 



e(n,/,T = 0) 
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The pressure of the baryons, at T = 0, defined as 

2 d(e/n) 



is given by 

P(n,I,T = 0) 
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In Eq. (JSOD E° F is the Fermi energy of symmetric nuclear matter corresponding to equilibrium 
density n and is given by 
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The Fermi momenta of protons and neutrons are 
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The chemical potentials of protons and neutrons, at T = 0, are given by 
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Figure 1: Energy per baryon of symmetric (a) and pure neutron matter (b) of the present model 
(MDIM) in comparison with those originated from realistic calculations. More details for the 
models UV14+TNI, UV14+UVII and AV14+UVII in Ref. [ffj and for the models A18+UIX and 
A18+du+UIX* in Ref. [22]. 
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Figure 2: a) The contribution of the various terms Va, Vb, Vc and the total potential energy 
density Vi n t as a function of the baryon density b) The momentum dependent term Vc as a 
function of the baryon density at temperature T = 0, T = 10 and T = 30 MeV. 



20 



40 




1 2 3 4 5 1 2 3 4 5 1 2 3 4 5 



k(fm" 1 ) k(fm') n(fm" 3 ) 



Figure 3: A comparison of the single-particle potential of symmetric nuclear matter from the 
present model (MDIM) with the microscopic calculations of Wirigna [37] and Li et at jH], for 
densities n = 0.1 fm -3 , n = 0.3 fm -3 and n = 0.5 fm -3 . 




Figure 4: The difference F(n, T, Y p ) — F(n, T,Y p = 1/2) as a function of (1 — 2Y P ) 2 at temperatures 
a) T = 10 and b) T = 30 MeV, for three baryon densities. 
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Figure 5: The Helmholtz free energy F(n, T, I) of /3-stable matter versus the baryon density n, 
for various values of T (in MeV) . 




Figure 6: (a) The free energy of symmetric nuclear matter (for T = 5 MeV and T = 20 MeV) and 
(b) for pure neutron matter (for T = 3 MeV and T = 20 MeV) of the proposed model (MDIM) 
in comparison with the free energy calculated by Friedman and Pandharipande model (FP) [7]. 
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Figure 7: The internal energy E(n,T) of /^-stable matter as a function of the baryon density n 
for various values of T. 




Figure 8: The internal energy (for T = MeV) for symmetric nuclear matter (SNM) and pure 
neutron matter (PNM) calculated with the MDIM in comparison with the FP model. 
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Figure 9: The nuclear symmetry energy calculated with the MDIM in comparison with the FP 
model as well as the results of Refs. [M], [58] and [55] . 
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Figure 10: The thermal energy E therrna[ (n,T) = E(n,T) — E(n,T = 0) of /3-stable matter versus 
the baryon density n, for various values of T. 
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Figure 11: The pressure P of asymmetric nuclear matter for Y p = 0.1 and 0.3 at a) T = 10 and 
b) 30 MeV. The full lines give the results calculated from Eq. (pHl) . while the squares represent 
results obtained by differentiating F(n,T) (Eq. ff23l ). 




Figure 12: The thermal pressure Pthermai( n , T) = P(n, T) — P(n, T = 0) of /3-stable matter versus 
the baryon density n, for various values of T. 
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Figure 13: The fractions of protons Y p electrons Y e and muons of /3-stable matter as functions 
of the baryon density n, for various values of T. 
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Figure 14: The Fermi-Dirac distribution function f T (n,T) for protons and neutrons (r = p,n 
respectively), for n = 0.2 fm™ 3 , n 



0.4 fm and n = 0.6 fm and various values of T. 
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Figure 15: The entropy per particle S of asymmetric nuclear matter with Y p = 0.2 at T = 10, 20, 30 
MeV. The full lines give the entropy calculated from Eq. (|2"TI) . while the squares give results 
obtained by differentiating F(n,T) (Eq. (1231) ). 
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Figure 16: Contributions to the total entropy per particle of protons (S p ) (up triangles) neutrons 
(S n ) (upside down triangles) and the total entropy (S^) (squares). 
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Figure 17: Contributions to the total entropy per particle of electrons (S e ) (up triangles) muons 
(Sfj) (upside down triangles) and the total Si (squares). 




Figure 18: The pressures of baryons Pb leptons Pi (electrons+muons) and the total pressure P 
versus the baryon density, n for various values of T. 
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